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Abstract 

A Landau-Ginzburg functional of two order parameters (charge-density (f) and 
mass-density deviation ij) is developed in order to yield a field theory relevant to 
ionic lattice g 3jSGS clS well as a family of off-lattice models of ionic fluids that go 
beyond the restricted primitive model (RPM). In a mean- field (MF) approxima- 
tion an instability of a uniform phase with respect to charge fluctuations with a 
wave-number k ^ is found. This second-order transition to a charge-ordered 
phase terminates at a tricritical point (tcp). Beyond MF, a singularity of a mass 
correlation function for A; — > occurs at ion concentration lower than that of the 
MF tcp. An effective functional depending only on r] is constructed. For low 
ion concentration the usual Landau form of the simple-fluid (Ising) functional is 
obtained; hence in this theory the critical point is in the Ising universality class. 

1 Introduction 

Over the past few years there has been growing interest in the critical behavior and 

phase separation of ionic fluids and of the Hamiltonian model that has been most widely 
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used in trying to understand the experimental results, the restricted primitive model 
(RPM), which we shall discuss below. The interest stems in part from the fact that on 
the experimental side there has emerged no clear picture of universality in the critical 
behavior of these fluids of the sort that describes a wide variety of non-ionic fluids. 
For example, in a series of impressive experiments done on a family of ionic fluids, 
Narayanan and PitzeiS have obtained results consistent with Ising-like behavior. But 
within the family they studied, they encountered a wide range of reduced temperatures 
t = {T — Tc)/Tc that mark the crossover from behavior well-described by mean-field 
theory to behavior that is Ising-like as one approaches the critical temperature T^. from 
above along the critical isochore. (Some of the crossovers come at much smaller t than 
encountered in the case of simple non-ionic fluids.) That means that not all members 
of the family can be well- described by the RPM; whatever the crossover temperature 
of the RPM, it is given by a single reduced temperature t according to the law of 
corresponding states. 

Theoretical studies have revealed no simple ionic-fluid universality either. For ex- 
ample the most comprehensive simulationi"@ and theoreticallHl studies both suggest 
that the critical point of the rpmHB is a simple liquid-gas critical point, whereas 
studies of lattice-gas versions of the rpmIiIB reveal that on simple cubic and body- 
centered cubic lattices one should expect tricriticality instead of simple criticality. This 
is in striking contrast to the situation one finds in non-ionic fluid models, where one has 
come to expect off-lattice and lattice models to share the same critical behavior. The 
shape of the phase diagram found in the continuous RPM from recent Monte Carlo 
simulations&fi is shown in Fig.l. The critical behavior revealed by these simulations 
is consistent with a critical point related to the density fluctuations that represents 
an instability with respect to fluctuations of wave number k = 0. In the case of the 
RPM with the positions of the ions restricted to the lattice sites, i.e., a lattice RPM, 
quite different phase behavior has been found, both in simulation and in mean-field re- 
sultsBlli as well as in a general theoretical study.0 From these references, what emerges 
is an instability of the uniform phase induced by charge fluctuations with wave num- 
ber k ^ 0. The line of the second-order transition to the charge-ordered phase with 
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the period 2T[/k just below the transition (A-hne) appears to terminate at a tricritical 
point, which hes in a different range of densities and temperatures than the critical 
point found in the continuum-space RPM studies, as shown in Fig. 2. 

In order to better understand this question of universality in ionic-systems we have 
developed a field-theoretic approach for studying the boundary of stability of a uniform 
phase in such systems. As prototypic microscopic models we consider both the lattice 
and the continuous versions of the RPM. In our treatment the differences between the 
lattice and the continuous systems are reflected in different forms of the couplings in 
the Landau- Ginzburg (LG) functional, which has the same general structure in both 
cases. With our method we are able to reproduce the essential properties of ionic 
systems both for the continuous-space and the lattice RPM within one formalism. We 
can also include the effects of the short-range interactions between various species, i.e. 
extend our model beyond the RPM. 

As in the earlier lattice RPM studies we have cited, an instability of the uniform 
phase with respect to charge fluctuations with 7^ is already found at the mean- 
field level when we apply our approach to the lattice RPM. We find second-order and 
first-order transitions to a charge-ordered phase at temperatures above and below that 
of a tricritical point respectively. The density at the tricritical point is found to be 
1/3 on the simplest mean-field level for both the lattice and the continuous versions of 
our LG functional. The structure of the charge-ordered phase depends on whether the 
positions of the ions are or are not restricted to the lattice sites. 

In addition to the instability induced by the charge-fluctuations with /c 7^ 0, an 
instability related to the mass-fluctuations with A; = occurs in a different region of 
thermodynamic space. This other instability, which has already been investigated on 
the mean-field level in the context of several standard ionic-fluid approaches applied 
directly to the continuum-space RPM&S can be detected when the charge-fluctuations 
are integrated out in our formalism. We find an equation for the critical point consistent 
with the fluid-theory result found previously by one of us.lll We also derive the effective 
functional depending only on density deviations from the average value, i.e. with 
the charge fluctuations integrated out. Near the critical point this functional has the 
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same form as the "0^" functional representing the Ising universahty-class. Hence in 
our treatment of the RPM model the critical point induced by the mass fluctuations 
belongs to the Ising universality class. 

Since we find the tricritical point observed previously only in the lattice RPM as well 
as the liquid-gas type critical point found previously only in the continuous RPM our 
approach enables us to give a unified description of the continuum-model and lattice- 
model behavior of ionic systems. Moreover, it sharply reveals two distict transitions 
that can be driven by the Coulombic forces in such models. One is a transition from a 
charge-disordered to a charge-ordered phase. The other is a liquid-gas type transition 
that does not involve charge ordering at all. In the continuum-space version one expects 
at least two solid phases associated with different charge ordering0 as well as the fluid 
phase, and it may be that the locus of charge-ordering transitions we find in our field 
theoretic versions of the RPM will occur in the continuum-space RPM as a solid-solid 
transition, or when one passes from a fluid to solid state. A reliable description of 
transitions involving solid phases of the continuum-space RPM that can be used to 
give a global-description of both its fluid and solid phases requires going beyond the 
treatment we have given here. 

In lattice-gas versions of the RPM, the details of the charge ordering that can occur 
will surely depend sensitively upon the lattice geometry and upon the extent of the 
core of the interparticle potential excluding multiple occupancy of neighboring sites 
as well as a central site as they do in non-ionic lattice-gas models. It seems likely 
that in some extended-core models, one will find both the liquid-gas-like criticality 
and phase separation characteristic of the continuum RPM as well as charge-ordering 
at higher densities involving phases that are also density ordered and hence solid-like. 
The lattice-gas version of the RPM and its extended-core generalization may well prove 
valuable in modelling the transitions one finds in the continuum-space RPM solid. 

Beyond the RPM we find evidence that the position of the tcp is extremely sen- 
sitive to the strength of the additional, short-range interactions, as for example the 
interactions between the ions and the particles of the solvent. For reasonable values of 
the parameter describing the effect of these interactions in a simple extension of the 
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RPM we find within our formalism that the tcp hes in the immediate vicinity of the cp 
found in the recent simulations, either slightly above or slightly below the cp in tem- 
perature, with its location depending very sensitively on the values of the parameter 
describing the effect of the short-range interactions. This raises the possibility that in 
some real ionic fluids in which there are substantial short-range interaction involving 
the solvent, one might flnd a tcp realizable in the liquid state or incipient as a tcp 
involving met ast able liquid phases. 

2 The functional in the case of RPM (no ion-solvent 
interactions) 

Our choice is designed to be relevant to a lattice-gas version of the RPM as well as 
variants thereof, such as the standard off-lattice RPM and an extension that includes 
the effect of ion-solvent interactions. 

In order to study criticality in ionic systems within a fleld-theoretic framework, 
we flrst construct an appropriate Landau- Ginzburg (LG) functional. We start the 
construction by calculating the grand thermodynamical potential Q^'^^ in the MP ap- 
proximation. For the mixture Q^^^ has the form 



In the lattice case the integration in (|2.1| ) should be replaced by a summation over the 
lattice sites, with the lattice constant a equal to the diameter of the ions, pa and pa 
are the density and the chemical potential of the component a respectively and the 
summation convention for the repeated Greek indices is used. JJ^^ is the energy in the 
MP approximation. is the Helmholtz free energy of the hard-core reference system, 
in the lattice (and ideal) case given by 




(2.1) 




(2.2) 



X 



In continuous models 




(2.3) 
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with 

f3j^ = -^-c\p). (2.4) 

OpaOpf^ pa 

Here we will limit ourselves to given by ( p.2| ). More accurate forms of in the 
continuous models will affect the explicit expressions for the coupling constants of the 
final LG functional, but not its general structure. 

For the RPM model a = +,—,0 correspond to a cation, an anion or a solvent 
particle (a 'hole') respectively, and Q^^^ can be expressed in terms of the local charge 
0(r) = p^—p^ and ion-mass p(r) = p^+p^ = 1—po densities. The energy contributions 
to Q^^^ read: 

in a continuum 



U^^ -fia [ dvpaiv) = 1 [dr [ dr%{r - r')0(r)0(r') -p [ drp(r) 

J ^ J J\r-r'\>a J 



(2.5) 



on a lattice 



U""^ - /i„ J]p«(x) = ^ E E ^^(^ - ^')0(x)0(x') - pEp(x), (2.6) 

X X x'^x X 

where /i = /i4. — /io = /i- — po- Vdr) = q'^/{Dr), with q being the charge and D the 
dielectric constant, is the Coulomb potential. Note the restrictions on the integration 
in (|2.5| ) and the summation in (p.6|) , which are necessary to avoid the self-energy 
contributions. These restrictions have a strong effect on the form of JJ'^^^ in the Fourier 
representation. It is convenient to introduce 

V{r) = Vc{r)e{r~a). (2.7) 



Then, with Vc replaced by V, the integration in (|2.5| ) and the summation in (|2.6| ) are 
over the whole space, with no restrictions, and 



jMF I d]i. ~. 



= J ^V(A:)0(k)0(-k). (2.8) 

The tilde refers to the Fourier transform of the corresponding function. In the contin- 
uum version, an upper cutoff ~ 27r/a is assumed as in standard theories, and in the 
lattice case the integral is strictly restricted to k = {ki, k2, k^) such that — vr < /cj < tt 
(fcj measured in units a~^). The explicit forms of V{k) are 
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in the continuum 



V{k) = in 



cos k 



(2.9) 



The restriction on integration in (p.5|) , |r— r'| > a, results in the factor cos k multiplying 
the Fourier transform of the Coulomb potential, 
on the lattice 



V{k) = 27T 



1 



3(1 -/(k)) 



with 



dh r dks 1 

° ' 27r-y_ 27r3(l-/(k)) 



1 

/(k) = -^cos/cj. 



(2.10) 



(2.11) 



(2.12) 



The first term in ( 2.10 ) is the Fourier transform of the Coulomb potential Vc and the 
second term results from the fact that the summation in ( |2.6| ) is restricted to x 7^ x'. 

The Q'^^^ assumes a minimum for the charge density = and the mass density 
p = Pq, and is next expanded in a power series in and 77 = p — pg. Assuming that 
and rj are both of small magnitude and vary slowly on the length scale of the lattice 
constant or particle core diameter a, we truncate the expansion and in the lattice case 
use a continuous approximation for fi*^^. 

The above procedure, described in more detail elsewhere,0 leads to a functional of 
the two-order parameter (OP) fields, and 77, of the following form: 



(2.13) 



with the Gaussian part given by 

f^^9 = ^J h'^'(r) + A2r]\r)] + 1 y dr y rfr'0(r)/3y(|r - r'|)0(rO (2.14) 



and with 



pa 



int 



dr 



(2.15) 
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In the above, r is a dimensionless distance, such that a = 1; h.o.t. stands for 'higher 
order terms'. In Ref. [17] we obtain for given by 



as = P*o-' (2.16) 



A2 = m-Pir' 



(2.17) 
(2.18) 



04 = 2p*-3 (2.19) 
where we use the standard notation in which pg is the dimensionless density pa^ and 

f3V{r) = f3*r-^e{r - 1) (2.20) 
with P* = T*~^ = I3q^ /Da. Higher-order couplings can also be expressed in terms of 

The functional defined in Eqs. ( |2.13|) -( p.20| ) is of lowest order in and 77, allowing for 
studying boundary of stability of the uniform phase in ionic fluids within bifurcation 
analysis. Without the h.o.t.., however, Vt becomes unstable for (f),ri —>■ 00. From 
Eqs ( |2.13| )-( |2?T5| ) one can easily see that the instability occurs if A = (63/8— 04^42)712 > 
[or explicitly pg < 1/3 from Eqs. (^TTB[)-(^rT^)] and corresponds to 0, 77 ^ 00 in such 
a way that r] = hcjP' with (-63/6 - \/~K)/A2 <b< (-63/6 + y/A)/A2. For one of the 
OP fixed, Q remains finite for the other OP going to infinity. When the h.o.t. up to 
0(77^) are included, Q becomes stable. Therefore at the later stage of our study (|^) we 
will also include the h.o.t. up to 0(?7^), which stabilize the functional. 

In the Fourier space Qg has the representation 

/3^, = IJ dk {c%ik)\m\' + Mmi') (2.21) 

where 

in the continuum 

Cl{k) = p;-'(l + An(3*pf-^^ (2.22) 
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on the lattice 



(2.23) 



As already discussed, in the lattice case the integral in (|2.21|) is restricted to — tt < 
ki < 71. 

In simple fluids the short-wavelength fluctuations are irrelevant in the critical region, 
and the long-wavelength fluctuations occur with the same probability in the continuous 
and in the lattice systems. This is basically the reason of the universality of critical 
phenomena. 

For k ^ i.e. for distances much larger than the lattice constant, both (|2.22| ) and 
( p.23| ) are of the same form: Ak^^ + B. For the long- wavelength charge fluctuations 
both the lattice and the continuous systems behave in a similar way, as in the case of 
the simple fluids. Note, however that the long - wavelength fluctuations are strongly 
suppressed. The Boltzmann factor ~ ea;p(— r/]), measuring the probability of 
fluctuations (0, approaches zero for charge fluctuations (f){k 0) [see Eqs. ( |2.21| )- 
( p.23|) ]. For a fluctuation 0(k) the charge neutrality is violated within regions of linear 
size 271 /k and the long- wavelength fluctuations are very rare. The relevant charge 
fluctuations are of flnite wavelengths. A probability of a fluctuation 0(k) with k flnite 
is different in the continuous and in the lattice case (compare ( |2.22|) and (|2.23|) ). In the 
RPM model the critical phenomena correspond to long-wavelength mass fluctuations, 
which are, however induced by the charge-fluctuations. Since the relevant charge- 
fluctuations are of short-wavelengths, and the probabilities of the short-wavelength 
fluctuations in continuum and on the lattice are different, one may expect differences 
between the lattice and the continuous systems. 



3 MF results 

Within the MF there is no ordinary critical point of a liquid-gas type, because the OPs 
are decoupled in the Gaussian part fig, and A2 never vanishes. Rather, an instability 
of the uniform phase with respect to charge fluctuations (pik) with k ^ occurs when 
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C^^{k) = 0. The boundary of stability of the uniform phase corresponds to the highest 
temperature for which C^^{k) vanishes. The bifurcation hues are 
in the continuum 

27rsinfcf, ^ 

T:{pI) = -^^Po- (3-1) 

where the wave number corresponding to bifurcation is given by 

2 

tankb = — — (3.2) 
kb 

from which one obtains kb ~ 2.46 and the slope of the bifurcation line ~ 1.61. 
on the lattice 

Tb* = 2n (Vo - ^) p; (3.3) 

and the wave vector at the bifurcation is 

kfe = (±7r,±7r,±7r). (3.4) 

The line of a continuous transition to a charge-ordered phase becomes first order 
at a tricritical point (tcp). To find the position of the tcp we first minimize Q with 
respect to r]. From the condition 6Q/6r] = 0, we can express rj in terms of 0, 

b 

r/(r) = r/o(r) = (3.5) 

and in the vicinity of the bifurcation we have 
1 r dk 



m 2 



(3.6) 

4 4 



4! J {27tYJ {2tiYJ {27tYJ {2ttY^ ' 



where 



A = a4 - J- = ?>pl-' [pl - , (3.7) 

e is the bifurcation parameter, e ~ ^/\^^'^^f\JT* , such that C^^{k) = 0{e'^). At local 
minima of Q { 6Q/6(f) = and 6Q/6ri = 0) = 0(e) and r] = O(e^). From the above 
we find that the tcp occurs at = 1/3 for the model given by Eqs. (|2.13|) - (|2.20|) . 
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Just below the bifurcation various metastable structures can be described as linear 
combinations of planar waves 

0(r) = $cos(kb ■ r) (3.8) 

where is a direction of oscillations and the amplitude is $ = 0(e). 
on the lattice 

r = {ri,r2,r^) are restricted to integer rj. For the integer rj (|3.8|) and (|3.4|) give 
(f){r) = ±$. Thus two sublattices, one positively- and the other one negatively charged 
occur. A charge-ordered phase with such a structure has been also found directly in 
the lattice RPM in MF results and in Monte Carlo simulations. 

in the continuum 

To determine the structure stable just below the continuous transition, one has to 
find for what linear combinations of planar waves (|3.8|) Q assumes the global minimum. 
For example, the structure of ionic crystal (like NaCl) corresponds to a superposition 
of four waves, with the four vectors k^ forming a tetrahedron.0 In Ref. [17] we verify 
by explicit calculations of Q approximated by ( p.7|) and (|3.7|) that for a single wave 
( p.8| ), i.e. for a lamellar structure, Q assumes a lower value than for the other struc- 
tures. Formation of charged layers, with the neighboring layers oppositely charged, is 
energetically favorable compared to the disordered phase when the average distance 
between alike ions within one layer is larger than the distance between oppositely 
charged layers. The ratio between the average distance between the ions within one 
layer and the distance between oppositely charged layers can be estimated on the basis 
of Po the period of modulations A = 2iT/kh ~ 2.55a. Near the tcp this ratio is 
> v/(2/A)3/p* ^ 1.2. 

We should stress here that for a — > (or with 6{r — 1) in (|2.2CI|) omitted), i.e. for 
point charges, no such instability takes place, because then C°^(fc) > 0, and consis- 
tently, A — > for a 0. Note also that the tcp lies far away from the critical point 
(cp), for which the latest simulation estimate iJl^ {p*,T*) ^ (0.08,0.05). 

The charge-charge correlation function in the MF approximation, G^^ = (5°^^ , 
never diverges for A; ^ 0, rather, G^^ — for A; — > (see ( [^.22| ) and ( ^^.23] )), i.e. 
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/ drG^^{r) = 0. By inspection of Eqs. (|2.22|) and ( |2.23| ) for ^ we can easily 



identify the inverse Debye screening length na = y^Anp^P*, describing the decay of the 
charge-charge correlations. 

4 Critical mass fluctuations 

Although there is no 'ordinary' criticality within the MF, the charge fluctuations may 
lead to the critical behavior of the mass-mass correlation function for A; — > 0, because 
the two OP-s are coupled beyond the Gaussian part of fi. Fluid theorieJUii and 
simulations&B predict low p* at critical; thus we concentrate on < p^^ = 1/3. 

In order to determine whether the criticality of the liquid-gas type occurs in this 
model, we calculate the correlation function for the mass fluctuations Gnrj, defined as 

G,,(ri,r2) = {v{r,Ur2)) - {viri)){v{r,)), (4.1) 

in standard perturbation expansion about the Gaussian solution. Next we examine the 
behavior of the Fourier transform of for — * 0, searching for singular behavior. 

The most probable mass fluctuation ?7o(r), accompanying a given charge fluctuation 
0(r), is given by (|3.5|) (if the h.o.t. in Eqs. ( p.l3|) -( ^T5D are included, rio{r) contains 
additional contributions of a form ~ 0^"; see Ref. [17]). Arbitrary mass fluctuation can 
be written as i] = r]Q + Ar] = —-^^(p'^ + Ar]. The functional i]] given by Eqs. ( |2.13| )- 
(|2.15|) can be expanded about the minimum at r/o- Then one obtains a functional of 
0, At] of the form 

f3n[<f),Ar]] = (3nauM + J dr{A7]{r)f + 0(Ar/^). (4.2) 
The fiaua;[0] IS glvcu by (^3), or in the real space representation by 

Pn,^, = jdv(^-j dv'Cl^{\v - r'|)0(r)0(rO + ^A0'(r) + 0(0*^)) (4.3) 



It is convenient to introduce the functional ( |4.2|) , because A2 is just a constant in this 
model and A?7(r) contribute only to local parts of the mass-mass correlation functions. 
The nonlocal correlations, relevant in the critical region, are determined by the auxiliary 



functional Vtaux ( |4.3|) , depending only on the single OP 
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In this description the one- and n-point correlation - functions for mass fluctuations 
can be expressed in terms of average values of the field 0^(r) and of a product of 0^ at 
n points, respectively. The explicit relations for the nonlocal parts of the correlation 
functions read 

(^(r)) = -^(0^(r))_ (4.4) 



2 



GZi^ur,)=[^] G^.^.(ri,r2) (4.5) 



where 

G^2^2{vu V2) = ((0'(ri) - (0'(ri)),„,)(02(r2) - (0'(r2))a«.))anx (4.6) 

and 

Gl[{r„...,r^)=(^-^^ G„^.(ri,...,r„) (4.7) 

In the above G'„02 is defined in a way similar to the definition of (j(^2^2 in Eq. ( [4.6|) . 
{...)aux means averaging with the Boltzmann factor ~ exp( 

In order to find the critical point we consider G^2^2, which is just proportional to 
G^lj, in the perturbation expansion, and verify whether there occur singular contri- 
butions for /c — i> 0. G'^2^2(ri, r2) can be expressed, in the standard way,0 in terms 
of Feynmann diagrams of the connected four-point function, with the vertex A4 and 
the free-propagator in which two pairs of external points are identified with each 
other. Next the Fourier transform for every diagram is calculated. Because G^^{k) is 
regular for k —>■ 0, and even tends to zero (see Eqs.( |2.2^ ) and ( |2.23D ), every individual 



diagram is regular for k —>■ 0; hence, any finite sum of diagrams is regular as well. 

From Eq. (|3.7|) we see that the effective coupling constant A4 is negative for 
Po < 1/3, i-e. for densities lower than the tcp density. The negative is of crucial 
importance for the behavior of G^2^2. For illustration consider a Fourier transform of 
a series of chains of loops, with the n-th order term (n vertices Aa and n+1 loops ) 
given by 

a^k) = i-A,nS^\k) (4.8) 
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with 

-goik) = \j drG%ir)V^^ = \ j dk'G°^(k')Gj/k - k'). (4.9) 
A contribution to ($0202 given by the sum of such diagrams has a form 

oo 

Y,^n{k)=Uk^)[l + AMk)r^- (4.10) 



In the part of {pq,T*) to which we restrict our attention here, ^o(^) is a smooth 
function (see ( ^ and (|2]2|), ( ^^231) ). For A < 0, however, (|]Tg) represents a 



singular contribution to G"' for A; if 



1+ A^o(O) = 0. (4.11) 

The above equation is the lowest order approximation for the critical singularity. A 
more accurate equation has a form 

1+A^(0) = 0, (4.12) 

where g{k) is a Fourier transform of the function representing a sum of all connected 
diagrams of four-point functions with two pairs of external points identified, which 
cannot be split into two distinct diagrams by splitting a single hyper-vertex ^4 into 
two two-point vertices. In other words, g{r) has no contributions which are of a form 
of chains; go is the lowest order approximation for g. For pQ ^ and with g{r) 
approximated by G^^{r)/2, where G,^^ is the charge-charge correlation function, ( |4.12| ) 
reduces to an equation obtained previously by one of us within fluid theory.0 

By summing up a related geometric series (as in Eqs.( [4.8D and ( [4.10| ), with go 
replaced by g) we obtain a more accurate form of the singular part of GJJ^, which we 
will denote by G^^, 

G',,{k) = m + A,mr' ■ (4.13) 

In the critical region we can expand G°^(/c) = G°^(/c)~^ about k = and we obtain 

C'^,ik)=ao + a2e (4.14) 
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with ao and a2 expressed in terms of g{0) and g{0)" in a standard way. 

Consider the higher-order mass - correlation functions and related vertex func- 



tions for n > 3 and with the two-point function approximated by ( [4.13|) 



G'::;(ri,...,r„) = -J dr'...J rfr"r„,(r',...,r") x ^^^(ri, rO...G;^(r„, r") (4.15) 



Unless the terms ~ 0^, ~ 0^ are explicitly included in ( p.l5|) , there are no contributions 



to Tn-q at the zero-loop level. At the one- loop approximation, for T^^, for example, we 
have 

r3,(ri, r2, rg) = G%{r,, r2)G'J^(r2, r3)G'5/r3, r^) ■ num.fac. + 0(2 - loop) (4.16) 
and similarly for n > 3. After integration ( |4.1(j| ) yields 



Tn = J rfri...t;r„r„^(ri,...,r„) ~ Gj/0)" = (4.17) 

Higher-order contributions to r„ also contain a factor (5,^^(0) = 0. Hence, for n > 3 
Tn give negligible contribution to the effective potentiahlli 

In the perturbation expansion for the charge-charge correlations G^^{k), every Feyn- 
mann diagram in the Fourier representation contains a factor G^{k) with n > 1, which 
for k = vanishes. In an infinite series of diagrams, similar to those which lead to 
divergent mass correlations, a segment of a corresponding chain for G^,^ is proportional 
to G^^{k). The infinite series of Feynmann diagrams gives a contribution to charge- 
charge correlations which vanishes for k ^ (compare ( |4.1CI|) with ^o(^) replaced by 
a term proportional to G^^{k)). Hence the charge-correlation length remains finite at 
the critical point, at which the mass-correlation length diverges. This is consistent 
with the results found in our earlier liquid-state theoretic analysis. 



5 Construction of the effective functional for < 

1/3 

In order to construct an appropriate effective functional depending only on rj in the 
critical region, we recall that Q given by ( |2.13| ) without the h.o.t. is unstable for (j),ri ^ 
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oo with ?7 ~ 0^ if Pq < 1/3 . We thus include the terms stabihzing Q and we consider 
the functional 



(5.1) 



with 



where, by construction of the functional 



(5.2) 



A, 



(5.3) 



A, 



1 1 ' 



if is approximated by 



The effective functional of r] is formally defined by 



Thus, 



(5.4) 



(5.5) 



(5.6) 



with 



-m<t>,v] 



(5.7) 



It is not possible to perform the functional integration in (|5.71) exactly, but it is 
also not necessary for the purpose of characterizing the mass-correlation singularity. 
We only need an approximate functional, which leads to the same behavior of the 
correlation functions at large distances in the critical region that the full JF does. The 
above-described analysis of the critical behavior of the correlation functions leads to 



dk 



Sri(k)5ri{-k) 



+ corr. 



(5.^ 
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where Crir^i^) given by ( |4.14| ) and 

Sf] = 7] — {rj) (5.9) 



with {rj) given by ( [4.41 ). With the above form of JFq the two-point mass correlation 
function ([4.13|) is recovered and (p) = p* + {i]) Terms ~ rj" with n > 3 are neglected in 
(p.8|). As usual, for small amplitude fluctuations such terms are negligible compared to 
terms ~ 7]'^, except for the critical fluctuations for which the second-order term vanishes 
at the critical point. For such fluctuations, i.e. ri{r) = const, however, the higher-order 
contributions to JFq vanish (unless the terms 0{(j)^),0{(f)^),0{(j)'^ri),0{(fP'7f) etc. are 
exphcitly included in r/], see (|4.17| )). 



By shifting and rescaling the field. 



i){v) = ( r/(r) + 4^ ) (5.10) 



we obtain from ( |5.6| ), ( |5.2| ) and ( |5.8| ) a familiar form of the Ising universality-class 
functional, namely 

a//M = j cir[-/.^(r) + ^^(r)2 + l(V^(r))2 + ^^(r)^] +/i.o.t. (5.11) 

with explicit expressions for the couplings corresponding to the functional ( |5.2| ) and 
(EH) given by 



^4 V 3^4/ 

U2 = (^ao - a^^ (5.13) 

and 

M4 = ^402'^ (5.14) 

Higher-order contributions to (|5.1| ) as well as more accurate approximation for the 
hard-core reference system (Eqs. ( |2.3| ) and ( [^.4D instead of (^^.2])) modify the couplings. 
In particular, to An {n = 3, 4) additional terms, related to the coupling constants of 
the terms ~ ^k^2n-k gj^ould be added, if the latter are included in the functional ( |5.1| ). 
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Such terms have significant effect on the position of the critical point. In the MF 
approximation for the functional ( p.ll| ) the cp is given by /i(pq,/3*) = U2{Pq, P*) = 0, 
which reduce to ^43 = 0:0 = 0. The density of the critical point, given by ^3 = 0, 
changes significantly when ^3 is replaced by a sum of several terms depending on 
p* resulting from the 'h.o.t.'. However, the 'h.o.t.' have no qualitative effect on the 
universality class. 

The U2 depends on P* through and 02, which in turn are functions of ^(0) 
and ^(0)" (see Eg. ( ^.14]) and below). The function g{k) is determined by the charge 
fluctuations, and for A; ^ depends on the whole spectrum of the fluctuations, as one 
can see from the lowest-order approximation for g, go, given by Eq.( ^4.9| ). Because of 
that dependence on the shape of G'^^{k), the position of the cp, U2 = 0, also depends 
on the form of G^^{k). In particular, for different systems, characterized in our theory 
by different forms of G^^ (as, for example, the continuous ( |2.22| ) and the lattice ( |2.23|) 
models) the relative locations of a possible cp signaling an Ising-like transition and a 
possible locus of transitions to a charge-ordered phase can be quite different. As a 
result, one or the other of these transitions can be preempted by the occurrence of the 
other. Thus the continuum-space RPM fluid appears to support only the cp while the 
cubic lattice-gas RPM appears to support only the charge-ordering transition with an 
attendant tcp. It will require analysis beyond that given here to locate the cp in our 
field-theoretic approach. 



6 Beyond the RPM (short-range interactions in- 
cluded) 

Here we limit ourselves to the continuous case. The lattice case is described in Ref. [17]. 
On our level of description the short-range interactions, including those between the 
ions and the solvent, can be represented by terms a20(V0)^/2 and A2^^(Vri)^ /2 as 
in standard Landau theories for short-range interactions. In general, 02 and A2 differ, 
respectively, from those of Eqs.( |2.l6| ) and ( ^.17| ), which refer to the absence of the 
short-range interactions. The and are the (MF) spin and mass correlation lengths. 
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respectively, in a system with the Coulomb interactions turned off. They should be thus 
comparable to the typical distance between the ions and the particles of the solvent. 
For our length unit both should be thus of order of unity. Here we restrict ourselves 
to the case of 02 and A2 given by (|2.16| ) and ( p.l7| ) and to a single length C, = C,^ = ^-q, 
to illustrate on the simplest level the crucial role of the interactions beyond the RPM 
on the position of the A-line and the tricritical point. 

The bifurcation analysis of the extended functional ( p.l3|) , which in addition con- 
tains the terms a2^^(V0)^/2 and A2^'^{'VriY /2, gives the generalization of Eqs. ( |3.1| ) 
and ( p. 21) for the bifurcation line 

and 

The tcp is very sensitive to the changes of ^ and is given by the analysis of Ref. [17] as 

* _ 1 

" 3 + sek! ■ ^^-^^ 

For ^ = 1, for example, = 1) = 0.026 and T;^p(^ = 1) = 0.007, which is well below 

the cp (compare the position of the tcp in the absence of the ion-solvent interactions, 
Plpi^ = 0) = 1/3 and T;,p(e = 0) = 0.54). 

In the case of ions large compared to the particles of the solvent, as for example in 
the system A^2226-B2226 studied by several groups^! a reasonable estimate for the typical 
distance between the ions and the particles of the solvent should be around 0.5. From 
( |6.3|) and Eqs. (|6.1| ) and (|6.2| ) we find the tcp for two slightly different values of ^, 
^ = 0.5 and ^ = 0.55 at (p*,T*) = (0.08,0.055) and (0.07,0.043) respectively. Note 
that in the first case the tcp lies slightly above the RPM cp=(0.08, 0.05), actually found 
in the simulations,&fl whereas in the second case the tcp lies below the observed cp. 

It is important to note the approximations and assumptions we are making in using 
Eqs. ( |6.1|) , (|6.2| ), and (|673|) as we do. There are precise values of and associated 
with each thermodynamic state in our system and, in particular, with the states along 
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the bifurcation line associated with ( |6.1| ), (|6^), and the tcp associated with ( |6.3|) . 
We do not know these precise values, however, so we make simple estimates of their 
magnitudes and then see how sensitive the tcp T* and p* that come out of (|6.1|), ( |6.2|) , 
and ( |6.3| ) are to slight changes in ^. We find great sensitivity, which suggests that the 
tcp location is very sensitive to the interaction parameters upon which the ^(f, and ^.^ 
depend. 

Whether a similar sensitivity of the cp location upon these parameter values exists is 
an important question that we cannot probe on the basis of our results here but which 
is deserving of further work. Such a sensitivity might help explain the remarkable 
disparities that different groupslll have found in the cp location in A'"2226-B2226 samples. 
An even more dramatic possibility is raised by the observation that comes out of our 
work that depending on the details of the short-range interactions in an ionic system 
one might have both a cp and a tcp [as in Fig. 3] or only one of these occuring in 
stable equilibrium, with the other describing metastable criticality that lies within the 
coexistence region defined by the stable singular point [as in Fig. 4] . We report further 
details in Ref. [17]. 
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Figure captions 



Fig. 1.- Liquid-gas coexistence curve of the continuous RPM model according to recent 
simulation studies.M Here Po is the dimensionless concentration of ions and T* the 
reduced temperature, defined after Eq. ( p.20| ). 

Fig. 2.- (a) A-line and coexistence curve from cell histograms taken from Monte Carlo 
simulation results for the lattice RPM desribed in Ref. [12]. (b) A-line and coexistence 
curve from the simple mean- field theory described in the Appendix of Ref. [12] for the 
same lattice RPM. (Note the difference in the temperature scales.) 

Fig. 3.-Generic phase-diagram associated with our field-theoretic description of the 
RPM, shown schematically. Here p and T denote the ion-concentration and tempera- 
ture respectively in arbitrary units. We show the case in which both a critical point 
and a tricritical point are realized, along with a A-line (dashed) representing a locus of 
transitions to a charge-ordered state lying below the A-line in temperature. 

Fig. 4.- Possible phase-diagrams associated with our field-theoretic description of the 
RPM with additional short-range interactions a2^^(V0)^/2 and A2^'^{VriY /2 included, 
shown schematically. Here p and T denote the ion-concentration and temperature 
respectively in arbitrary units. We show the two special cases described in the text: 
(a) the critical point is realized only as a metastable point, (b) the tricritical point is 
realized only as a metastable point. 
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Ciach-Stell, Fig. 3 
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Ciach-Stell, Fig. 4(a) 
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